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Final state interactions at the threshold of Higgs boson pair production 
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We study the effect of final state interactions at the threshold of Higgs boson pair production in 
the Glashow-Weinberg-Salam model. We consider three major processes of the pair production in 
the model: lepton pair annihilation, ZZ fusion, and WW fusion. We find that the corrections caused 
by the effect for these processes are markedly different. According to our results, the effect can cause 
non-negligible corrections to the cross sections for lepton pair annihilation and small corrections for 
ZZ fusion, and this effect is negligible for WW fusion. 


Half a century after the construction of the Glashow- 
Weinberg-Salam (GWS) model of electroweak interac¬ 
tions, the observation of Higgs boson at the Large Hadron 
Collider (LHC) is its latest triumph [ 2 , @. In the future, 
further determination of the properties of Higgs boson 
will be extremely important for us to understand the fun¬ 
damental laws of nature. Any deviation from the predic¬ 
tions in the Standard Model can give us a valuable clue 
to the long-hunted new physics. In order to ascertain the 
role of Higgs field in the Standard Model, we need to mea¬ 
sure its couplings to fermions and gauge bosons. Further¬ 
more, to reconstruct the details of the Higgs potential, 
we have to precisely determine the Higgs self-couplings. 

At present, an important task of us is to explore poten¬ 
tial effects for the Higgs self-interactions, and the Brout- 
Englert-Higgs mechanism must be precisely tested from 
low to high energies. In this paper, we shall investigate 
the effect of final state interactions at the threshold of 
Higgs boson pair production in the GWS model. This 
effect is tied to the low-energy properties of the Higgs 
self-interactions. 

After the electroweak symmetry breaking, the Higgs 
self-interactions can be written in the form 

V(H) = \m\H 2 + XvH 3 + ±H 4 . (1) 

For the mass of Higgs boson m# ~ 125.5 GeV and the 
vacuum expectation value v ~ 246 GeV, A = m 2 H /2v 2 ~ 
0.130. 

To establish the non-relativistic potential for the Higgs 
self-interactions, we need to consider the scattering pro¬ 
cess HH —>• HH to leading order in A. There are four 
Feynman diagrams that contribute, see Fig. |I] 

Note that we can separately discuss the contributions 
from each diagram, since there is Bose statistics involved 
in the process. In the beginning, let us consider the con¬ 
tributions from the quadrilinear self-coupling of Higgs 
boson. The diagram is shown in Fig. [lji, and its ampli¬ 
tude is 


iM = —*6A. 


( 2 ) 



FIG. 1. The tree-level diagrams for the scattering of two Higgs 
bosons in the Standard Model. Solid line denotes Higgs boson. 


According to the definition of Born scattering ampli¬ 
tude in non-relativistic quantum mechanics, the non- 
relativistic potential in the momentum spaces is 

F quad (q) = 2^T- ( 3 ) 

Using the Fourier transform to V(q), we can get the non- 
relativistic potential for this term 

3 A 

FjuadM = ^r<5 (3 V). (4) 

The presence of the delta function comes from the spe¬ 
cially local structure in the diagram. 

The amplitude for the “annihilation” diagram in 
Fig. [TJd is 


iM = — 18 \itih 
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where p\ and P 2 are the 4-momenta of the incoming Higgs 
bosons. The amplitude in the non-relativistic domain can 
be approximated as Eq. ©, and then we can get a well 
defined non-relativistic potential. It is not necessary to 
repeat the same calculations. The non-relativistic poten¬ 
tial for the “annihilation” diagram is defined as 
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The amplitude for the “exchange” diagram in Fig. [lj; 
is 


iM = — l8\m 2 H -- - 5 -, (7) 

q z — m H 

where 4-momentum q is the momentum transfer. In the 
non-relativistic domain, the amplitude becomes 

iM = 18 A m 2 H -—- 2 ”, ( 8 ) 

I q?+m 2 H 

and the non-relativistic potential in momentum space is 


V ex (q) = 


9A 1 
2 | q | 2 +m 2 H ' 


(9) 


After the Fourier transform, the non-relativistic potential 
for the “exchange” term is 


Kx(r) = (10) 

87rr 

The presence of the “cross” term in Fig.QJi comes from 
the indistinguishable property of the identical particles. 
However, the “cross” term for the scattering of two iden¬ 
tical particles does not contribute to the non-relativistic 
potential. It is because that the principle of identity 
also needs to be considered in the scattering theory of 
non-relativistic quantum mechanics, and then the non- 
relativistic potential that we get from the elastic scat¬ 
tering of the distinguishable particles will automatically 
include the contributions from the “cross” term of the 
identical particles. In appendix, we shall show an ele¬ 
mentary example in detail. 

As a result, we obtain the non-relativistic potential for 
the Higgs self-interactions 

V(r) =Fquad(r) + Vann ) + V ex {r) 

=^-5 {3 \r)~-e- mHr , ( 11 ) 

m H r 

where coupling constant a = 9X/(Sn). 

We know that the probability of scattering is in pro¬ 
portion to the squared modulus of the wave function of 
the created particles and a non-relativistic potential in 
the final state can alter the amplitude of the final state 
wave function in the reaction zone Q- Notice that Higgs 
boson can be treated as a stable particle for the Higgs 
self-interactions, since the width of Higgs boson is only a 
few MeV pj. Therefore, the non-relativistic Higgs poten¬ 
tial in the final state might be able to modify the relevant 
cross section in perturbation theory, and this effect can 
be very important for us to understand the low-energy 
properties of the Higgs self-interactions. After getting 
potential m, we shall investigate the effect of final state 
interactions at the threshold of Higgs boson pair produc¬ 
tion. In recent years a similar effect, which considering 
the distortion of the incoming wave function of the initial 
particles, in dark matter was discussed, see. e.g., IMS- 



FIG. 2. Higgs boson pair production in the processes of lepton 
pair annihilation (a-c), and the generic diagram of the Higgs 
self-interactions in the final state (d). 


At first let us consider the final state interactions in 
the processes of lepton pair annihilation // —>■ HH, see 
Fig. [2] We know that in perturbation theory the scatter¬ 
ing amplitude is defined in the plane wave bases. How¬ 
ever, as we mentioned above, due to the appearance of 
the interactions in the non-relativistic region of the pro¬ 
cesses, the wave functions in the reaction zone can be 
very different from the plane waves. To consider the cor¬ 
rections in the processes, the universal amplitudes for the 
s-, t-, and M-channel diagrams in Fig. fl2J can be written 
in the forms 

Ms = M*M0)r f (0)/$m?(0), (12) 

Mt = Mt'tpi(rt)'ipf(rt)/'ipi{r t )'tp ( f* (r t ), (13) 

M u = MW(a)^(a)M 0 (a)i(/*(a), (14) 

where t/>» and ipf are the wave functions for the relative 
motions of the particles in the initial and final systems, 
ip® = e lki ' r and ip® = e lkf ' r are the plane waves for 

free motions in the initial and final systems, r t and r u 

denote the separate distances of the two vertexes in t - 
and u-channel diagrams, and M p denotes the relevant 
amplitude in perturbation theory. 

The wave functions introduced in universal amplitudes 
(11211141) can give the corrections to the squared modulus 
of the plane waves in the reaction zone, and we can solve 
the Schrodinger equations to get the wave functions if we 
know the forces between the particles in the initial and 
final states. Notice that in the perturbation region the 
universal amplitudes become the same as the amplitudes 
in perturbation theory. 

In the process of the electron-positron annihilation, it 
is straightforward to see that the contributions of t- and 
u-channel diagrams in Fig. ([2]) can be ignored, since the 
Yukawa coupling to electron ( m e /v ) is very weaker than 
the Higgs self-couplings. Therefore, in practical appli¬ 
cations, the squared modulus of the amplitude at the 
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threshold can be simplified as 


\M\% 


ee^HH 




(15) 


where |'i/ ; }(0)| 2 is the correction term caused by the final 
state interactions in the process. 

Note that to find the corrections caused by the ef¬ 
fect of final state interactions, we only need to consider 
the corrections to the s-wave, because the Higgs boson 
pair is emitted in an s-wave at the threshold. The re¬ 
duced two-body radical Schrodinger equation with a po¬ 
tential V(r) is d 2 cj)(r)/dr 2 — mV{r)(j>{r ) = — (mu) 2 (/>(r), 
where <f>(r) = rRo(r ), and Rq is the s-wave radical func¬ 
tion. The radical wave function R^i can be normalized 
as / 0 °° r 2 Rk'iRkidr = 2ir6(k' — k), and notice that Ru 
is real in this convention. Since we have already ob¬ 
tained non-relativistic potential (fill) , to find |^>i(0)| 2 in 
Eq. m we can numerically solve the Schrodinger equa¬ 
tion with boundary conditions —> 0 as r —> 0 and 

<j>(r) —> 2sin (fcr + Sq) as r -> oo. However, we should 
note here that although the non-relativistic potential for 
the Higgs self-interactions is composed of different terms, 
only the “exchange” term would be involved in the effect. 
The reason is that the final state interactions between 
two separate particles come from the exchange of their 
force carriers, thus only the “exchange” part of the non- 
relativistic potential need be considering. The numeric 
simulations for the 125.5 GeV Higgs boson are presented 
in Tab. Q] We find that at the threshold the two-Higgs- 
boson final state interactions can increase the cross sec¬ 
tion by near ten percent. 

TABLE I. The numeric simulations for the corrections in 
electron-positron annihilation process, where y/s is the center 
of mass energy. 


y/s{ GeV) 

Correction(%) 

251.1 

9.84 

251.2 

9.83 

251.3 

9.82 


The strength of the effect grows slowly close to the 
threshold, and it is interesting to investigate that how the 
largest corrections can be reached. Using the principle of 
detailed balancing Q 


1 da 


Pf 


dQ 


f) = 


- —(f 

p 2 dn [J 
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where the states i* and f* are the time-reversed relative 
to the states i and /, and the momenta p a = m a v a (a = 
i,f.) for the reduced masses and the relative velocities of 
the two-body systems, we can construct the connection 
between the cross sections for the diagrams in Fig. [2] and 
its inverse process 




a 


HH-yff 


Is — 4 m 2 H 
8 s — 4 to 2 


(17) 


Note that here we have already considered the spin statis¬ 
tic weights and the symmetry factors for the processes. 
This relation holds true at all energies. 

The r.h.s of Eq. ( 1171 ) would become the square of the 
Higgs boson velocity vjj at non-relativistic energies, and 
we notice that the general theory of scattering requires 
that the cross section for the process HH —> ff is di¬ 
rectly proportional to v^ 1 in the low-energy limit and the 
other component of the cross section would be a constant 
0 . Using Eq. (El), we can find that a should be 
in direct proportion to vh in the low-energy limit. 

To the leading order in perturbation theory the cross 
section = \ 2 m 2 VH/ (167 tto|j) at the threshold. 

We here do not show the full expression which contains 
the extremely small contributions of the t- and u-channel 
diagrams, because it is cumbersome and no more instruc¬ 
tive. It can be seen that the cross section has already got 
a factor vh in the expression. Consequently, we can con¬ 
clude that the corrections from the diagram in Fig. [2ji 
must be independent of the kinetic energies in the low- 
energy limit. The numeric simulations in Tab. [TT] show 
that the largest corrections to the cross sections cannot 
exceed ten percent. 


TABLE II. The corrections for the “static” Higgs boson pair 
production in the processes of electron-positron annihilation. 


^(GeV) 

Correction(%) 

251.001 

9.85 


It should be noted that, as a matter of fact, a routine 
tree-level calculation can show that at the threshold s- 
channel diagram always absolutely dominates in the lep¬ 
ton pair annihilation processes. Hence, the results given 
above can also apply to muon and tau leptons. 






FIG. 3. Production of a Higgs boson pair in WW/ZZ fusion 
processes. 

The processes of WW/ZZ fusion are important for the 
Higgs boson pair production in the GWS model, and 
now let us consider the final state interactions in these 
processes. The diagrams for the tree-level processes are 
shown in Fig. [3] 
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The universal amplitude for the processes is 

M = A4q U ad + A i s + Mt + Mu, (18) 

where A4 qua d is the universal amplitude for the diagram 
in Fig. [Si. 

Notice that the contributions from the t- and it-channel 
diagrams can not be ignored in the processes. Thus, to 
find the corrections caused by the final state interactions, 
besides ipf( 0 ), we have to find ipf( r t) and V’/( r «) in the 
Mt and Mt- This fact causes a serious problem for us 
to calculate the corrections caused by the effect of final 
state interactions. In principle, we cannot obtain the 
exact numeric simulations for the corrections, because 
it is impossible to get the r t and r u in the diagrams. 
However, there is a standard approach in quantum theory 
that can help us to offer the first theoretical predictions 
for the effect in WW/ZZ fusion processes. Notice that 
the energy of the non-relativistic Higgs boson E ~ mu 
and the uncertainty principle in the relativistic case is 
Ar ~ 1/E [TTj | . It indicates that in the diagrams the 
distances around l/rriH are important for the system of 
the created low-energy Higgs boson pair. To calculate 
the effect of final state interactions, we can assume that 
the r t and r u in the diagrams are 1 /m//, and then the 
numeric simulations could approximate the strength of 
the corrections at the threshold. 

Before doing the calculations, we should note that in 
the lepton pair annihilation processes we calculate the 
corrections of the squared modulus of the wave function 
at the origin, and then without any specific considera¬ 
tion there is only the s-wave of the ip® that naturally 
contributes. However, here we shall not calculate the 
corrections to the s-wave at the origin in the Mt and 
M u . Thus we need to resolve the plane wave as e lkfZ = 
S“ 0 (— i) 1 {2l + \)Pi{cos 9 ){t / k f) 1 {d/rdr) 1 sin kfr/kjr and 
then take account of the s-wave alone. 

In the calculations we take Mz = 91.188 GeV and 
Mw = 80.385 GeV fl[, and we calculate the effect close 
to the threshold. The numeric simulations for the cor¬ 
rections in ZZ fusion process are presented in Tab. m 
We find that the strength of the effect for ZZ fusion is 
one order of magnitude weaker than that for lepton pair 
annihilation. The corrections in this process are small, 
but the effect might need to be considered for the high- 
precision experiments. 


TABLE III. The numeric simulations for the corrections in 
ZZ fusion process. 


Vs(GeV) 

Correction(%) 

251.01 

0.853 

251.001 

0.855 


The numeric simulations for WW fusion are shown in 
Tab. IIV1 We find that the strength of the effect for WW 
fusion is even one order of magnitude weaker than that 
for ZZ fusion. In general the effect in this process is 


negligible. Meanwhile, we notice that the numeric sim¬ 
ulations for WW fusion are sensitive to the distances of 
the interactions in the t- and it-channel diagrams, and in 
the future further studies on the accurate theoretical pre¬ 
dictions for the corrections in WW fusion process could 
be interesting. 


TABLE IV. The numeric simulations for the corrections in 
WW fusion process. 


Vs(GeV) 

Correction(%) 

251.01 

0.082 

251.001 

0.084 


In summary, we have calculated the effect of final state 
interactions in three major processes of Higgs boson pair 
production for the first time, and we find that the correc¬ 
tions caused by the effect for these processes are markedly 
different. This paper could be useful for understanding 
and precisely determining the Higgs self-interactions at 
low energies. 
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APPENDIX 

The “cross” term for the identical particles does not 
contribute to the non-relativistic potential. In this ap¬ 
pendix, we shall give an elementary example in detail. 

Let us consider the scattering of two electrons in QED, 
see Fig. [I] 


(a) (b) 

FIG. 4. The scattering of two electrons in QED. 

The cross section (M/ller formula) is 

da a 2 (2 E 2 — m 2 ) 2 

dH qed 4 E 2 {E 2 — m 2 ) 2 

4 3 {E 2 — to 2 ) 2 4 

X W0 " ihTd + (2E 2 — m 2 ) 2 1 + 

(19) 

where m is the mass of electron. 
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In the non-relativistic domain, the cross section be¬ 
comes 


It is elementary to find that the cross section for the 
process is 


7 NR 2 i i 1 

da a 1 1 1 

dflQED 16m 1 2 3 u 4 5 sin 4 |0 + cos 4 \6 sin 2 A0cos 2 i0 

( 20 ) 


where v is the velocity of the electron. 

It is well-known that the Coulomb potential a/r comes 
from the diagram in Fig. [3Jr for the distinguishable 
charged particles. Now let us consider the scattering 
of the identical electrons via the Coulomb potential a/r 
in non-relativistic quantum mechanics. The principle of 
identity requires that the asymptotic wave function must 
be 


v> = e ikz ± e~ ikz + -e ikr [f(6) ± /(tt - 0)]. (21) 

r 


da 

dQ QM 


a 


16 m 2 v 4 sin 4 


cos 4 \6 


1 Oi 1 

—-7— cos( —— In tan 2 -0)1. 

sin 2 i0cos 2 ^0 y 2v 2 u 

( 22 ) 


In the perturbation region of the scattering process, 
Eq. (l22l) becomes the same as Eq. © , see e.g., 

The example given above illustrates that the non- 
relativistic potential that we get from the “distinguish¬ 
able” term will automatically take account for the contri¬ 
butions from the “cross” term for the identical particles. 
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